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Surface Bound States and Spontaneous Current in Cyclic rf-Wave Superconductors 
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On the basis of Eilenberger theory, surface bound states and spontaneous current are studied in cyclic rf-wave su- 
perconductors as a broken time-reversal symmetry of superconductivity in cubic lattice symmetry. We discuss how the 
spontaneous current and the electronic states depend on the orientation of the surface relative to the symmetry of super- 
conductivity. The condition for topological Fermi arcs of zero-energy surface bound states to appear is identified in the 
complex pairing fimction of the cycUc d-wave. 
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Among unconventional superconductivity, the broken time 
reversal symmetry (BTRS) of superconductivity is an im- 
portant topic of study, because it has exotic properties, such 
as zero-energy surface bound states and spontaneous mag- 
netic field by spontaneous current. In most of the previous 
studies of the BTRS superconductor, the pairing symmetry 
was assumed to be a chiral p-wave px + ipy, considering the 
cases of the superfluid ^He A-phase'"^ or Sr2Ru04.'''^ The 
spontaneous magnetic field of a BTRS superconductor is de- 
tected by muon spin rotation (yuSR) experiment in Sr2Ru04,'* 
PrOs4Sbi2,*'^ LaNiC2,** PrPt4Gei2,'' and LaNiGa2."' In a yuSR 
experiment, PrOs4Sbi2, PrPt4Gei2, and LaNiGa2 show dif- 
ferent types of relaxation curves from those of Sr2Ru04 and 
LaNiC2. Therefore, there may be some variety in the types 
of BTRS. Since the crystal lattice symmetry of PrOs4Sbi2 is 
cubic Oh (or more exactly, Th), the pairing symmetry may be 
different from the chiral /7-wave of Sr2Ru04. Therefore, it is 
important that we examine the possibility of a new type of 
pairing symmetry of BTRS other than the chiral />-wave, and 
study the properties of the new pairing function to identify the 
pairing symmetry. 

In this study, we consider the BTRS of superconductivity in 
cubic lattice symmetry. From the classification table of possi- 
ble pairing symmetries under spin-orbit coupling on the basis 
of point group theory," as a BTRS state keeping Oh sym- 
metric superconductivity, we find -H cjj-ky + <jj|fe? in spin- 
singlet pairing and kxX H- (jj±kyy H- <xi\k^z in spin-triplet pair- 
ing with = e±i2;r/3 ^jjg former is cyclic cZ-wave pairing, 
and it is stable in weak coupling theory. The latter is non- 
unitary spin-triplet pairing, and it is not stable in weak cou- 
pling theory. Thus, the case of cyclic J-wave pairing is studied 
here. The possibility of cyclic cZ-wave superconductivity and 
its exotic properties have been discussed in the study of Bose- 
Einstein condensation^^'* and fermionic superfluidity'^ in 
cold atomic gases. There, the non-Abelian ^-fractional vor- 
tex is also possible in an isotropic atomic gas system. Thus, 
the BTRS of cyclic J-wave pairing is one of the interesting 
pairing symmetries for the theoretical study of topological su- 
perconductivity. 

Studies of the surface state are important to identify the 
pairing symmetry of unconventional superconductivity. For 
example, in a <i;c2_y2-wave superconductor, the surface bound 
state depends on the surface orientation, and zero-energy sur- 



face bound states appear for the (1,1,0) surface."^' The char- 
acteristic dispersion relation of the surface bound state reflects 
the pairing symmetry of the bulk superconductivity. In the 
dxy-^3Nt pairing and the /5 j-wave pairing, we have flat disper- 
sion of the surface bound states for the (1,0,0) surface. '^'^'' In 
superfluid ^He, surface bound states show topological Fermi 
arcs for the A-phase and a Majorana cone for the B phase.^'^ 
Therefore, studies of the surface bound state are also neces- 
sary for the complex pairing function of cychc c?-wave super- 
conductivity. 

The purpose of this study is to clarify the properties of sur- 
face bound states in a cyclic J-wave superconductor, as an- 
other example of the BTRS state, on the basis of Eilenberger 
theory. Since the structure of surface bound states depends 
on the relative angle of the surface and the pairing symmetry, 
we study the surface-orientation dependence of spontaneous 
current and electronic states in the surface bound state. We 
also discuss the condition under which the zero-energy sur- 
face bound state appears in the cycUc <i-wave superconduc- 
tivity. 

In our study, wave numbers are denoted as (ka,kb,kc) for 
crystal coordinates. Under an Oh crystal field, we assume the 
pairing function to be the cyclic J- wave given by 

0±(k) = a/I (^c + ^±kl + ojlkl) 

= ^(^(3^-l)±i^(^-^)). (1) 

Therefore, 0-(k) = 0+(k). In Eq. (1), k is mapped on the 
Fermi sphere and normaUzed aslc^-i-k^-\-k^ = 1. From Eq . ( 1 ) , 
the cyclic cZ-wave is a combination of c/j.2_^,2-wave and c/3j2_i- 
wave components. The amplitude |^±(k)| has Oh symmetry 
with 8 point nodes in (1,1,1) and equivalent directions. The 
phase of 0±(k) has Th symmetry. For 120° rotation around the 
(l,l,l)-axis (ka -> kb, h -> kc, kc -> ka), <p+ -» o/(/)+. For 
90° rotation around the (0,0,l)-axis (ka — > kb, kb — > -ka), 
(p+ — > <p-(i^ (p+). In the calculation of the surface state, we 
use the coordinates (kx, ky, k^) and (x, y, z), where the x-axis is 
fixed to be perpendicular to the surface at x = 0. {kx,ky,k^ 
is obtained by rotational transformation from (ka, kh, kc). The 
x-direction is defined as (cos cos 0, cos sin 0, sin 0) in the 
crystal coordinates. 
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Surface states are calculated on the basis of Eilenberger 
theory, following the method used to study the surface states 
of superfluid ^He.^ Since we consider the case of spin-singlet 
pairing here, the transport-like Eilenberger equation is re- 
duced to the 2x2 matrix form 



-iv ■ Vg = 



ico„ 
A*(k, r) 



-A(k,r) 

-iw„ 



for quasi-classical Green's functions 



^(k, r, w„) i/(k, r, w„) 
-i/"'"(k,r,dL»„) -g(k,r,u„) 



(2) 



(3) 



with the normalization condition = -n^i. cj„ is the Mat- 
subara frequency. We consider the case where an external 
magnetic field is not applied. For simplicity, we neglect the 
contribution of vector potentials. The pair potential takes the 
general form 



A(k,r) = A+(r)<^+(k) + A_(r)0-(k). 



(4) 



We treat A+ as the dominant component and A_ as the in- 
duced component near the surface, so that A_ = and A+ 
is a bulk value when x — > oo. Since we assume a spheri- 
cal Fermi surface, the normalized Fermi velocity is v = k. 
We assume that the surface condition is specular, and a quasi- 
particle with (-fc.v> K) is reflected to {k^, ky, k.) at x = 0. We 
solve the Riccati equation derived from the Eilenberger equa- 
tion [Eq. (2)], and obtain the quasi-classical Green's func- 
tions. Energy, temperature, and length are in units of Trk^Tc, 
Tc, and = Hvp/lnksTc, respectively. The pair potential is 
calculated by the gap equation 



Ai(r) = io2r£<<^:(k)/), 



(5) 



ai„>0 



where <■ ■ ■ )k indicates the Fermi surface average. Aq - Nogo 
is the dimensionless pairing interaction defined by the cutoff 
energy oj^ as 1 //Iq = InT + 2T YjZ^o '^n ■ ^® carry out cal- 
culations using the cutoff = ^^k^iTf^. The calculations of 
Eqs. (2) and (5) are iterated at T = ^.TT^ until self-consistent 
results are obtained. 

The current density of spontaneous current is given as 



J(r) = ;o2r2<vImg)k, 

M,,>0 



(6) 



and produces the spontaneous magnetic field, jo - SineNo/c, 
and A^o is the density of states at the Fermi level. The total 
current of the surface bound states is given by the integral 



(7) 



Jtotai = I J(r)dx. 
Jo 

The local density of states (LDOS) is obtained as 

N(E, r) = {N(E, k, r))u = No{Rsg(k, r, iw„ + i/7))k, (8) 

where we use the solution of the Eilenberger equation [Eq. 
(2)] for the real energy E in the self-consistently obtained pair 
potential. We use a small smearing of // = 0.007. N(E, k, r) is 
the k-resolved LDOS. 

First, in Figs. 1(a) and 1(b), we present the self-consistent 
solution of the pair potential and spontaneous current for 
some surface orientations. At the surface region, when the 
length is of the order of the coherence length, A+ is sup- 
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Fig. 1. (Color online) (a) Depth dependence of pair potential |A+(x)|, 
|A_(a')|, and spontaneou.s current \J(x)\ as a function of distance x from the 
surface. Solid lines are for the surface orientation (0,6) = (45°, 0) and 
dashed lines are for ((j), 8) = (20°, 0). |A±| is in the unit of nk^Tc- \3{x)\ and 
Uiotall are, respectively, normalized by the values |J(j: = 0)| and Utotall for 
(0,e) = (45°, 0). (b) Same as (a), but {0,9) = (45°, 25°) for solid lines and 
(0, 9) = (45°, 45°) for dashed lines, (c) Surface orientation (0, 9) dependence 
of the surface state |A+(x = 0)|, |A_(,r = 0)|, |J(;c = 0)| (solid lines), and 
IJiotiill (dashed lines). The horizontal axis indicates that the surface orienta- 
tion changes from (1,0,0) to (1,1,0) [6> = and = ^ 45°] and from 
(1,1,0) to (0,0,1) [0 = 45° and = ^ 90°]. (d) Schematic plot of spon- 
taneous current at the surface of a large spherical sample. The line indicates 
the path of the horizontal axis in (c). The path (1,1,1)-(0,0,1) is equivalent to 
the path (l,l,l)-( 1,0,0). 
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Fig. 2. (Color online) LDOS N{E, r) at surface x = with various surface 
orientations, (0,0) = (0,0), (25°, 0), (45°, 0), (45°, 25°), and (45°, 45°). 



pressed, and the induced A_ accompanies the y-component 
of the spontaneous current. These behaviors depend on the 
surface orientation. For the (1,1,0) surface with (0,0) - 
(45°, 0), since the t/v--y--wave component in the pairing func- 
tion changes sign in the quasi-particle reflection at the surface, 
|A+| = |A_| at the surface so that the c/j.2_,,2-wave component 
vanishes there. On the other hand, for the (1,0,0) surface with 
{(f>, 9) - (0, 0), A_ = and J = even near the surface, and 
the surface state is the same as the bulk state. 

To confirm how the surface state depends on the surface ori- 
entation, in Fig. 1(c), we plot the surface values \AJ,x - 0)|, 
|J(jc = 0)1, and Utotall along flie trajectory (1,0,0)-(1,1,0)- 
(1,1,1)-(0,0,1). Along (1,0,0)-(1,1,0), |A+(x = 0)| decreases 
and |A_(.Ji: - 0)| increases. |J(x = 0)| increases and is maxi- 
mum at = 20°. However, the total current Utotall monotoni- 
cally increases until = 45°. This is because the surface re- 
gion in which spontaneous current appears becomes wider as 
(p approaches 45°, as seen in Fig. 1(a). Along (1,1,0)-(1,1,1), 
|A-t(x - 0)1 and current decrease. For the (1,1,1) surface of 
the point node direction, A_(x = 0) = and 7(x) vanishes. 
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Fig. 3. (Color online) k-resolved LDOS N(E,k,r) at surface jc = as functions of ky and E for L = 0, 1/(2 V3), 1/ V3, and 0.75. k, = (1 - /c^ - kl)"^. 
Various surface orientations, (4>, 0) = (0, 0) (a), (25°, 0) (b), (45°, 0) (c), (45°, 25°) (d), and (45°, 45°) (e), are presented. 



Along (1,1,1)-(0,0,1), which is equivalent to (l,l,l)-( 1,0,0), 
|A+(x = 0)1 increases monotonically. |A_(x - 0)| and current 
decrease after they increase. 

In order to understand the flow of spontaneous current 
J, we consider the surface current in a large spherical su- 
perconductor of cyclic c/-wave pairing. As in Fig. 1(d), at 
the point node direction (1,1,1) and the equivalent points, 
J = 0. Around the node points (1,1,1), (1,-1,-1), (-1, 1,-1), 
and (-1,-1, 1), J flows clockwise. Around other node points 
(-1,-1,-1), (-1,1,1), (1,-1,1), and (1,1,-1), J flows 
counterclockwise. The direction of the flow is related to the 
angular momentum of the Cooper pairs via the phase winding 
of the spherical harmonic function, F2,m, on the Fermi sphere. 
The pairing function is expressed as (p+ oc 

F2,-2 - i V2F2,i 

when z \\ (1, 1, 1) and as oc ^2,2 + V2}'2,-i when z \\ 
(1, -1, 1).'^ The spontaneous current at the surface is induced 
by the change in the angular momentum of Cooper pairs when 
they are reflected at the surface. By summing the sponta- 
neous currents around neighboring point nodes, J is enhanced 
at (1,1,0) and the equivalent positions. On the other hand, 
at (1,0,0) and the equivalent positions, J is canceled to zero 
by summing neighboring spontaneous currents. These current 
distributions may be reflected in the spontaneous magnetic 
field distribution that is expected to be observed. 

Next, we discuss electronic states at the surface. In Fig. 2, 
we present the LDOS N{E, x = 0) for some cases of surface 
orientation. The LDOS is expected to be observed by tunnel- 
ing spectroscopy at the surface. For the (1,0,0) surface with 
{(p, 6) - (0, 0), surface states have the same electronic states 
as those in the bulk. There, A^(£') oc at low E owing to 
point node excitations. The gap edge at E ~ 0.8 comes from 
the maximum of |A+0+(k)| at (ka, kb,kc) oc (1, 0, 0). The hump 
at £■ ~ 0.4 corresponds to the saddle points of |A+(^+(k)| at 
(k„,kb,k,)^ (1,1,0). 

With increasing <p from (1,0,0) to (1,1,0), low-energy sur- 
face bound states appear, including zero-energy states. The 



peak at £■ ~ 0.8 decays, and a new peak appears at £ ~ 0.5. 
In the cases of surface orientations {4>,6) = (45°, 25°) and 
(45°, 45°), the LDOS N(E,x = 0) has high intensity at low 
energy, because both A+{x = 0) and A-(x = 0) are largely 
suppressed at the surface. 

In order to understand the structures of the surface bound 
state and spontaneous current, we study the surface orienta- 
tion dependence of the k-resolved LDOS N{E, k, r) at the sur- 
face X = 0, as presented in Fig. 3, to examine the dispersion 
relation of surface bound states. For the (1,0,0) surface with 
(0, 9) - (0, 0), the surface state is the same as the bulk state. 
Thus, N{E, k, X - 0) shows the gap structure of the pairing 
function |A+0+(k)|. In the panel for k~ - 0, we see the maxi- 
mum of the gap |A+0+(k)| to be £ ~ 0.8 at A:,, = 0, ±1, and the 
saddle point energy E ~ 0.4 at ky = + 1 / V2. In the panel for 
k. = 1/ V3, zero energy states appear due to the point nodes 
at A:,, = ±1/ V3. 

With increasing from (1,0,0) to (1,1,0), in-gap states of 
the surface bound states appear. In the panel for k. - 1/ V3, 
flat dispersions of zero-energy modes exist at 0.28 < \ky\ < 
0.77 for the surface with (0,0) = (25°, 0). For the (1,1,0)- 
surface with (0,0) = (45°, 0), zero-energy modes extend to 
all ky at k, - 1/ V3. This is a property of the iij.2_,,2-wave 
component, because the d^^i^i-wave component vanishes at 
k^ - \l Vs. In Figs. 3(b) and 3(c), the flat dispersions are 
raised to finite energies at k^ + 1 / V3. When k- > 1 / V3 (k- < 
1 / V3), they appear at £ > (£ < 0) for A;^, > and at £ < 
(E > 0) for ky < 0. Since a negative E indicates an occupied 
state, the imbalance of negative E states of ky > and ky < 
induces the spontaneous current. Since the contribution of the 
negative E state of ky > at L < 1 / V3 is stronger than that 
of ky < at L > 1 / V3, this diff'erence causes the spontaneous 
current to flow in the positive y direction. 

When the surface orientation changes from (1,1,0) to 
(0,0,1) upon increasing 6 at <p = 45°, the flat dispersions 
in Fig. 3(c) become dispersive, as shown in Figs. 3(d) and 
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Fig. 4. (Color online) Bold lines indicate wave numbers {ky, k,) of topo- 
logical Fermi arcs satisfying the condition for zero-energy states, Eq. (9). 
ky^ = (\ — k^, — i?)''^. Various surface orientation cases, (0,0) = (0,0), 
(25°, 0), (45°, 0), (45°, 25°), (45°, 34°), (45°, arcsin(l/ V3)), (45°, 36°), and 
(45°, 45°), are presented. The thin lines k. = 0, 1/(2V3), 1/V3, and 0.75 
correspond to the horizontal axis in the panels in Fig. 3. The intersections of 
bold and thin lines indicate that the zero-energy states appear at {ky, k-) in 
Fig. 3. 



diagonal lines. Zero-energy states also appear on the verti- 
cal line ky - 0. The endpoints of the zero-energy states on 
the vertical line are the point node directions. When ap- 
proaches arcsin(l/V3) ~ 35° of the (l,l,l)-surface direc- 
tion, all vertical and diagonal lines pass through the center 
{ky, kz) = (0, 0). The reconnection of these lines occurs at this 
point-node wave number, as shown in Fig. 4. From the panel 
for (0, 9) - (45°, 25°), we find that zero-energy states appear 
at ky = in the plots along the Hnes fc. = 1/(2 V3), 1/ V3, 
and 0.75. From the panel for (0, 9) - (45°, 45°), zero-energy 
states are seen to also appear at ky in the plots along 

= and 1/(2 V3). These well explain the wave numbers of 
the zero-energy states in Figs. 3(d) and 3(e). 

In summary, we studied the surface orientation dependence 
of the surface states in the cyclic d-wave superconductor, 
as an example of the BTRS state in cubic lattice symmetric 
superconductivity. There, spontaneous currents flow around 
each point node direction if we prepare spherical samples. We 
also identified the condition under which topological Fermi 
arcs of zero-energy surface bound states appear in cyclic d- 
wave superconductors. These are useful results for studying 
new types of BTRS superconductors other than those with 
chiral p-wave pairing. 

We thank T. Mizushima and T. Kawakami for fruitful dis- 
cussions. This work was supported by KAKENHI Grants No. 
21340103 and No. 24840048. 



3(e). We see the reconnection of the dispersion curve between 
= and 25° in panels for k, = 1 /(2 V3) and 0.75. The recon- 
nection also occurs between 9 - 25° and 45° in the panel for 
k, - 0. Zero-energy states appear at ky = in the upper three 
panels in both Figs. 3(d) and 3(e), indicating flat dispersion on 
part of the line ky = 0. In Fig. 3(e), we also see zero-energy 
states at ky in the panels for k, - and kr - 1/(2 V3). 

Lastly, we discuss the condition of k for zero-energy sur- 
face bound states to appear in the complex pairing function 
of cyclic (/-wave superconductivity. We define the phase cpi^ 
of the pairing function before reflection at the surface as 
(p+{-kjc, ky, k.) = \(l)+(-kx, ky, A;j)|e"^'° and the phase (/?ref after 
reflection as (f>+(ky,ky,kj - \(f>+(kjc,ky,k^)\&'^'" . As the condi- 
tion under which the zero-energy surface bound states in Fig. 
3 are well explained, we find the relation 

i.e., the ;7r-phase shift i^ref - fin = ^ (mod 2;7r).'^ Equation (9) 
can be applied even when {(p+i-k^-, ky, k~)\ + |</i+(A:v, ky, A:,)|. 

Wave numbers k satisfying Eq. (9) are presented in Fig. 
4. These topological Fermi arcs are terminated at the point 
node directions.^ For the ( 1 ,0,0)-surface with (</>, 0) - (0,0), 
there are no zero-energy surface states. For the (l,l,0)-surface 
with ((/>, 0) - (45°, 0), all ky values satisfy Eq. (9) on the Unes 
ky = ±1/ V3. When changing ^ from (1,0,0) to (1,1,0), the 
region of ky in which zero-energy surface bound states appear 
increases. These reproduce the behavior of the zero-energy 
flat dispersion in the panels for k, - 1/ V3 in Figs. 3(a)-3(c). 
When changing 9 from (1,1,0) to (0,0,1), the lines of zero 
energy at ^- = 1/ V3 are shifted to smaller and become 
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